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Abstract 

A new model for the dynamical structure of galaxies is studied based on the so 
called Liquid-Droplet model describing the structure of the nucleus. It is argued that 
in the galactic scale the inverse squared law of long range gravitational interaction may 
be replaced by a new law of gravitational interaction which is effectively short range 
while representing the characteristic features of Mach principle. Application of this new 
interaction at galactic scale gives the rotation curve in good agreement with observations. 
Then, the Virial theorem for this kind of gravitational interaction is developed and also 
the Tully-Fisher relation is obtained. Finally, a physical explanation is given for the 
so called constant acceleration in the MOND as the effective gravitational strength of 
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galaxies. 
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1 Introduction 

It is well known that classical Newtonian dynamics fails on galactic scales. There is astronom- 
ical and cosmological evidence for a discrepancy between the dynamically measured mass-to- 
light ratio of any system and the minimum mass-to-light ratios that are compatible with our 
understanding of stars, of galaxies, of groups and clusters of galaxies, and of superclusters. 
Observations on the rotation curves have turn out that galaxies are not rotating in the same 
manner as the Solar System. If the orbits of the stars are governed solely by gravitational 
force, it was expected that stars at the outer edge of the disc would have a much lower orbital 
velocity than those near the middle. In fact, by the Virial theorem the total kinetic energy 
should be half the total gravitational binding energy of the galaxies. Experimentally, however, 
the total kinetic energy is found to be much greater than predicted by the Virial theorem. 
Galactic rotation curves, which illustrate the velocity of rotation versus the distance from the 
galactic center, cannot be explained by only the visible matter. This suggests that either 
a large portion of the mass of galaxies was contained in the relatively dark galactic halo or 
Newtonian dynamics does not apply universally. 

The dark matter proposal is mostly referred to Zwicky [1] who gave the first empirical 
evidence for the existence of the unknown type of matter that takes part in the galactic scale 
only by its gravitational action. He found that the motion of the galaxies of the clusters 
induced by the gravitational field of the cluster can only be explained by the assumption of 
dark matter in addition to the matter of the sum of the observed galaxies. Later, It was 
demonstrated that dark matter is not only an exotic property of clusters but can also be found 
in single galaxies to explain their flat rotation curves. 
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The second proposal results in the modified Newtonian dynamics (MOND), proposed by 
Milgrom, based on a modification of Newton's second law of motion [2]. This well known 
law states that an object of mass m subject to a force F undergoes an acceleration a by the 
simple equation F = ma. However, it has never been verified for extremely small accelerations 
which are happening at the scale of galaxies. The modification proposed by Milgrom was the 
following 

F = mfi( — )a, (1) 
a 



1 if x > 1 

(2) 

x if <C 1, 



where a = 1.2 x 10 _10 ms -2 is a proposed new constant. The acceleration a is usually much 
greater than a for all physical effects in everyday life, therefore /i(a/a )=l and F = ma as 
usual. However, at the galactic scale outside the central bulge where a ~ a we have the 
modified dynamics F = m(^). Using this new law of dynamics for the gravitational force 
F = we obtain 

GM ,a\ 

— = - , (3) 

r l a 

which results in the constant rotational velocity 



v 2 = ^fGMa . (4) 



Another interesting model in this direction has been recently proposed by Sanders. In this 
model, it is assumed that gravitational attraction force becomes more like 1/r beyond some 
galactic scale[3]. A test particle at a distance r from a large mass M is subject to the acceler- 
ation 

a =^T9(r/r ), (5) 



where G is the Newtonian constant, r is of the order of the sizes of galaxies and g(r/r ) is a 
function with the asymptotic behavior 



Dark matter as the manifestation of Mach principle has also been considered as one of 
the solutions for the dark matter problem. According to Mach principle [4] the distant mass 
distribution of the universe has been considered as being responsible for generating the local 
inertial properties of the close material bodies. Borzeszkowski and Treder have shown that 
the dark matter problem may be solved by a theory of Einstein-Mayer type [5]. The field 
equations of this gravitational theory contain hidden matter terms, where the existence of 
hidden matter is inferred solely from its gravitational effects. In the nonrelativistic mechanical 
approximation, the field equations provide an inertia-free mechanics where the inertial mass of 
a body is induced by the gravitational action of the cosmic masses. From the Newtonian point 
of view, this mechanics shows that the effective gravitational mass of astrophysical objects 
depends on r such that one expects the existence of new type of matter, the so called dark 
matter. In what follows we will introduce a new model for the dynamical structure of a galaxy, 
based on a model for the structure of the nucleus, as an alternative to dark matter and MOND 
to explain the rotation curve. 

2 Liquid-Droplet model for the nucleus 

In nuclear physics, Weizsacker's formula or the semi-empirical mass formula (SEMF) is used 
to approximate the mass and various features of atomic nucleus [6]. The theoretical basis of 





1 if r > r 



(6) 



5 



this formula is the liquid drop model. This model in nuclear physics, first proposed by George 
Gamow and developed by Niels Bohr and John Archibald Wheeler, treats the nucleus as a 
drop of incompressible nuclear fluid which is made of nucleons (protons and neutrons) binding 
together by the strong nuclear force. It predicts the binding energy of a nucleus in terms of the 
numbers of protons and neutrons it contains. This statement manifests in the semi-empirical 
mass formula which has five terms on its right hand side. These correspond to i) volume binding 
energy of all the nucleons inside the nucleus, ii) surface binding energy of all the nucleons on 
the surface of nucleus, iii) the electrostatic mutual repulsion of the protons, iv) an asymmetry 
term and v) a pairing term (both derivable from quantum mechanical considerations). 

The first two terms concern about the attractive strong nuclear force. We shall focus on 
these two terms and leave the other terms irrelevant for our purposes. Assuming an approx- 
imately constant density for the nucleus one can calculate the nuclear radius by using that 
density as if the nucleus were a drop of a uniform liquid. The liquid droplet model of the 
nucleus takes into account the fact that the forces on the nucleons on the surface is different 
from those associated with the interior ones. This is something similar to taking into account 
of surface tension as a contributor to the energy of a tiny liquid drop. 

In the case of volume binding energy, when an assembly of nucleons of the same size 
is packed together into the smallest volume, each interior nucleon has a certain number of 
other nucleons in close contact with it. Actually, the number of pairs that can be taken 
from A particles is A(A — l)/2, so one may expect a term for binding energy proportional 
to A 2 interacting pairs. However, the strong force has a very limited range, and a given 
nucleon may only interact strongly with its nearest neighbors. In other words, the nuclear 



6 



force is saturated and each nucleoli contributes an almost constant energy to the binding of 
the nucleus. Therefore, the number of pairs of particles that effectively interact is roughly 
proportional to A and not A 2 , and the nuclear binding energy is proportional to the volume 
of the nucleus. The volume term suggests that each nucleon interacts with a constant number 
of surrounding nucleons, independent of A. 

In the case of surface binding energy however, each nucleon at the surface of a nucleus 
interacts with fewer other (interior and surface) nucleons than one in the interior of the nucleus. 
Therefore, the binding energy associated with each surface nucleon is less than that of the 
interior one. In fact, the surface term makes correction to the volume term by its negative 
contribution as follows 



where B accounts for the binding energy and a v , as are some appropriate constants with the 
same order of magnitude. One may focus more on the volume term to make maximum use of 
its physical justification. The physics behind this term in the case of a real liquid-droplet is 
as follows: To evaporate one liquid-droplet we have to use a definite heat energy as Q v M rn A 
where Q v accounts for the vaporization latent heat, M m is the mass of each molecule and A 
is the total number of molecules in the droplet. This heat energy is necessary to overcome all 
the attractive interactions between the molecules. Hence, it is exactly equal to the binding 
energy of the droplet. This justifies the first term as 



which results in the identification ay = Q v M m . Therefore, we come to the conclusion that 
the binding energy per each molecule, namely (B/A), corresponding to the volume term is 
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B = a v A - a s A 2/s , 



(7) 



By = Q v M m A = ay A 



(8) 



independent of the total number of molecules A. This has a simple reason: the number 
of nearest neighbours of each interior molecule is independent of the total size of the system 
(droplet), then (B/A) becomes independent of A. This is characteristic feature of every system 
in which the range of interaction between the constituent particles is short in comparison with 
the size of the system. Of course, this feature undergoes a surface correction due to the 
finite size of the system according to (7). Moreover, if we consider the repulsive coulomb 
force between the protons, then its negative contribution proportional to Z 2 (Z is the atomic 
number) affects the ratio (B/A) in a more complicate way. One may add also the contributions 
of an asymmetry term and a pairing term. In general, all of these contributions to the volume 
term except for the pairing term are negative. 

The result is that the curve of (B/A) with respect to A is almost flat except for the small and 
large values of A. The decrease for small values of A, is due to the fact that the full binding of 
each nucleon is accomplished when that nucleon is fully encircled by other neighbor nucleons. 
However, this does not apply for the surface nucleons; then for light nucleus having larger 
portion of surface nucleons the binding energy for each nucleon in the nucleus is decreased. 
On the other hand, the decrease for large values of A results from the coulomb repulsive force 
between each protons which is proportional to Z 2 and hence increases more rapidly than A, 
resulting in a decrease of binding energy per each nucleon. It is common to associate the 
almost flat portion of the curve with the saturation property of the nuclear force which itself 
is emerged due to the short range of this force. If this force would not be saturated and the 
nuclear force would be a long range force, then each nucleon would interact with all other 
nucleons with a term proportional to the size of the nucleus or A 2 . The more large size of 
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the nucleus, the more stabilized system, a result which is not favored by the nature. Hence, 
it seems the saturation property or short range of the nuclear force is related to the fact that 
nature avoids the structures with unlimited large size. 

3 Liquid-Droplet model for a typical galaxy 

To model the structure of the galaxy on the basis of Liquid-Droplet model we need to revisit 
the concept of the range of gravitational interaction in a large size structure of a galaxy. In fact, 
both the range of an interaction and the size of structure in which such interaction is governed 
are deeply related to the causal structure of space-time. Newton believed in the action of 
distance so that every gravitational effect, whatsoever far from a point, can be felt immediately 
at that point. The concept of a local gravitational field obeying the Lorentz invariance however 
changed this view: no gravitational effect at a large distance can be experienced at a point 
unless after a time interval that a signal with light velocity takes to enter the future light cone 
at that point. Even, in the Machian viewpoint where it is assumed that each gravitational 
constituent feels gravitationally all other constituents, the concept of local gravitational field 
obeying the Lorentz invariance is considered. This lorentz invariant causal structure then pose 
a crucial limit on the effective range of gravitational interaction in a large size structure. In 
fact, it is reasonable to assume that the effective range of any interaction in an structure is 
modeled by the size of that structure. For example, the nucleus with its small size determines 
an effective short range of Yukawa nuclear force. A deeper analysis on the basis of causal 
space-time structure may reveal that the large size of a galaxy determines the effective range 
for gravitational interaction inside this structure. Suppose, in a similar way that one calculates 
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the binding energy in a nucleus by Liquid-Droplet model, we calculate the binding energy in 
a galaxy. In the case of nucleus, we learned that each nucleon has full interaction just with 
neighbouring nucleons, due to short range of nuclear interaction. In the case of galaxy, we 
know gravitational interaction has long range. But, if we vaporize a small portion of a galaxy 
to free all the stellar mass in this portion from their mutual gravitational interactions, then 
all other stellar mass far enough from this portion of the galaxy will not feel this vaporization 
until its signal (with light velocity) can be reached by them. This exactly means that each local 
portion of a galaxy has an effective short range interaction compared with the large size of the 
whole galaxy. In other words, each small portion of the galaxy has just effective short range 
gravitational interaction with its limited neighborhood portions due to the limited velocity 
of light. This feature in the galaxy is similar to the feature in a massive nucleus where the 
range of effective nuclear force between close nucleons is short in comparison to the size of 
the nucleus. Therefore, if we accept in the context of Liquid-Droplet model that the effective 
range of gravitational interaction in a galactic structure is modeled by the effective size and 
mass of that structure, (an example of which is the short range nuclear force in the nucleus), 
then we may propose a model of gravitational interaction for the galaxy as follows 



where U is the effective gravitational potential in a galaxy containing N gravitational con- 
stituents of average mass m, R is the radius (effective size) of the galactic bulge , G is the 
gravitational constant, and k(r) is the r-dependent number of neighbor constituents for each 
individual constituent interacting with them. 

This form of gravitational potential may be partly interpreted as Machian . We know the 




(9) 
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mass is a positive quantity and it is not possible to screen the gravitational effects. There- 
fore, we assume that in any two (close) body process in the galaxy the extra influence of all 
other matter distributed over the effective size of the galaxy should be taken into account. 
Accordingly, we can replace the whole gravitational effect of all effective interactions at the 
position of a constituent of mass m by the effect of a spherical shell of the effective mass M 
and the effective radius R. This shell acts like a gravitational Faraday cage inside of which a 
gravitational potential exists 

. = (10) 

From the Liquid-Droplet model, we shall assume that each constituent gravitationally interacts 
with some limited neighbor constituents k. But, this number will depend on the distance r 
of that constituents from the center of the galaxy, as we shall explain later. By this kind of 
interaction we mean a local interaction. On the other hand, we assume a Machian (or nonlocal) 
interaction, as explained above, contributing to $ along with the local interaction. Therefore, 
if we consider the contribution of both local and nonlocal interactions to the gravitational 
potential (10) as the manifestation of effective mass M, and compare it with the gravitational 
potential energy (9) trough U = m$, then we obtain the effective mass M as an active 
gravitational mass 

M(r) = k(r)Nm. (11) 

This is a key feature of the present model in which the active gravitational mass is different 
from the additive mass of the galaxy, namely Nm. Moreover, we may associate the mean 
contribution —Gm 2 /R to each pair of constituents so that the potential energy at the location 
r will contain k(r)N terms where k(r) and N account for the local and nonlocal sectors of 
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interactions, respectively. 

We expect that the linear dependence of U on N will result in the saturation of the 
gravitational force, similar to the saturation of the nuclear force in the nucleus. In fact, 
similar to the volume binding energy of a nucleus which is calculated using the Liquid-Droplet 
model, the volume binding energy of the galaxy is given 



where the vaporization latent heat ay per each constituent in the galactic bulge is given by 



4 Rotation curve 

In this section, we study the rotation curve of a typical galaxy upon the Liquid-Droplet model. 
For each gravitational constituent located at a distance r from the center of galaxy having 
equilibrium state, one may approximately relate the order of magnitudes of the kinetic and 
potential energies as 



Now, we introduce the very important role of the variable k in the potential energy U. In the 
Liquid-Droplet model, k is the number of neighbor constituents interacting effectively with 
the individual constituent in consideration. In the case of a real Liquid droplet or a nucleus 
this number is a constant which is determined by the geometric shape of packed molecules or 
nucleons. However, this is not the case for the galaxy where there are large voids between the 
constituents. Assuming a spherical symmetry for the high density mass distribution in the 



B v = a v N, 



(12) 




(13) 




(14) 
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central galactic bulge, it is easy to realize that the number of neighbor constituents around 
each constituent i located in a volume element grows with the distance r similar to the volume 
element dV = 4nr 2 dr. Therefore, in the central galactic bulge we have 

k in = ar 2 , (15) 

where rj is the distance of the i th constituent from the center of galaxy and a is a constant. 
The potential (9) is then rewritten as follows 

U(n) = -Arl (16) 

where A = aGNm 2 /R is another constant. By substituting (16) into the equation (13) we 
obtain 

l -mv1 ~ Arl (17) 

which leads to the familiar result in that section of rotation curve which corresponds to the 
mass distribution in the galactic bulge, as 

Vi ~ n. (18) 

For the outside of the central bulge where the mass is smoothly distributed with very low 
density (compared with the bulge) over each volume element, it is plausible to assume that 
each constituent located in the middle of a volume element is surrounded by some of the other 
interacting constituents whose number is a very slow function of the distance r\. A suitable 
choice is 

k out = \nr u (19) 

for which we obtain 

U(n) = -B\nn, (20) 
13 



^mv?~B]nri, (21) 

where B = GNm 2 /R, or 

Vi ~ ( 22 ) 

This result gives an almost flat rotation curve for the outside of the central bulge, in complete 
agreement with current observations. 

We may repeat the calculations in a more dynamical way to obtain the almost same results 
by resorting to Newton's second law. In the case of galactic bulge the Newton's law becomes 

mvf dU 



n dr-i 

which gives the velocity proportional to the distance as 



2Ar t , (23) 




2A 

Vi = \—n. (24) 



m 



For the outside of the central bulge the Newton's law becomes 



mvf dU _ B 



n dr-i n ' 

which gives a constant velocity 



(25) 



v l = \-. (26) 



m 



In order to match the velocities at = R we need to equate the followings 



m V m 



which results in 



a = W (28) 
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Substituting this into Eq.(15) leads to the important result that each constituent inside the 
central bulge is at most (r; = R) effectively interacting with one half of another neighbor 
constituent. This shows that in the Liquid-Droplet model each constituent is locally interacting 
with a fraction k of other neighbor individual constituent, while globally interacts with TV 
distant constituent. This fraction is very small close to the center of galaxy but approaches to 
"one half at the radius R of the bulge. 

Note that all the above results come out due to our special and simple assumptions about 
the forms of the quantity k(r). It is worth noting that the quantity k(r) inside the galactic 
bulge seems to be almost the same for all galaxies, but outside the galactic bulge it may depend 
on some characteristic features of each galaxy like its shape, size, ellipticity, mass distribution 
, etc. Therefore, it is expected that in this model the rotation curve for the outside galactic 
bulge is not universally flat for all the galaxies, rather it differs for different galaxies according 
to their characteristic features; a result which is observed in some galaxies. 

5 Virial theorem 

For a system of N point particles, the virial G satisfies the following equation [7] 



where T is the total kinetic energy of the system, Fj is the net force on the i th particle at the 
position i-j, and G is defined as 



dG 
~dt 



N 



2T + ^F i -r,. 



(29) 



i=i 



N 



(30) 



i=l 
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We have realized that in the Liquid-Droplet model, the total force on each particle % is acted 
by a fraction of only a neighbor particle. Moreover, this fraction depends on the distance of 
that particle from the center of galaxy. Therefore, the mutual interaction between all particles 
which is usually considered in a gravitational system has now been replaced by this new type 
of limited interaction. If we assume the force to be derived from a potential energy, then 
we should accept that this potential is just a function of the distance of each particle from 
the center of galaxy r\ rather than the relative distance between the particles, namely 
Fj = —V r JJ. Therefore, we have 

N N dU 

EF^ = -£^r, (31) 

i=l i=l % 



Hence, Eq.(29) becomes 



dG _ "dU 

— = 2T-Y —n. 32 

dt dn y ' 



For the central galactic bulge with U in ~ rf we obtain 

2T-2J2 U tn (r t ) = 2T - 2U in , (33) 
dt ^ 

where Ui n is the total potential energy of the system in the galactic bulge consisting of N in 
particles. On the other hand, outside the central bulge of the galaxy consisting of N out particles 
with U out = —Blnri, we have 

dG Nout 

= 2T+J2B = 2T + N out B. (34) 
at i=1 

The virial theorem states that the time average (dG/dt) vanishes, hence we obtain 

(T) m = (U T ) in , (35) 
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for the inside, and 

(T)out = ~N out B, (36) 

for the outside of the galactic bulge . Therefore, the total time average of the kinetic energy 
of the whole galaxy is obtained as 

(T)tot = (U T ) m - l -N out B (37) 
GNm 2 fe ,, , 1 



R W 2 
By interpreting N out B as {U^out, we obtain 

(E) in = (T) m + (U T ) in = 2(T) m , (38) 

(E) out = (T) out + {U T )out = -(T) ou t, (39) 

(E) tot = 2(T) m -(T) out . (40) 

In a stable gravitationally bound system we require (E) tot < which leads to 

{T) m < l -{T) mt . (41) 

This is an interesting result which states that the total average kinetic energy in the galactic 
bulge is less than the half total average kinetic energy in the outside of galactic bulge. 

6 Tully-Fisher relation 

The Tully-Fisher relation is an empirical relationship between the intrinsic luminosity (propor- 
tional to the stellar mass) of a spiral galaxy and its amplitude of its rotation curve (which sets 
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the total gravitational mass) [8]. It has been turned out that on large scales most astronomical 
systems like a spiral galaxy have much larger mass-to-light ratios than the central parts. This 
fact may be explained by the relation (11) as follows. In the present model, the gravitational 
mass is the active mass M. On the other hand, the luminosity which is proportional to the 
stellar mass becomes proportional to the additive mass Nm. Therefore, the mass-to-light ratio 
becomes proportional to 

M 

k(r). (42) 



Nm 

Since in both the galactic bulge and outer region, k(r) is an increasing function of the distance 
from the center of galaxy r (see Eqs.(15), (19)), hence we have much larger mass-to-light ratios 
than the central parts which is in complete agreement with the Tully-Fisher relation. 

From another point of view we may study the Tully-Fisher relation in the present model. 
We may rewrite the equation (26) as 



m 2 - 



If we consider R as the effective size of the galaxy and Nm as the approximate stellar mass, 
then assuming a roughly constant mass density we may write Nm ~ R 3 or R 2 according 
to whether we assume a spherical shape or disk shape for the whole galaxy, respectively. 
Therefore, Eq.(43) reads 

v a ~ Nm, (44) 

where a = 3 or 4 corresponding to spherical or disk shapes for the galaxy, respectively. This 
is an interesting result which again accounts for the empirical Tully-Fisher relation. 
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7 Liquid-Droplet model and MOND 

According to MOND, the rotational velocity is given by v — (GNma )^. However we find 
that in the present model this velocity is obtained by v = ( °^ m ) 2 • Comparison of these two 
velocities leads to the following result 

GNm 

a = (45) 

This means the constant acceleration a which has no clear physical explanation in the MOND 
is explained in the present model through (45) indicating that this acceleration is nothing but 
the effective gravitational strength of each galaxy. To justify this relation we consider for 
example the Andromeda galaxy whose stellar mass Nm and bulge radius R are approximately 
7 x 10 U SM and 10 5 LF, respectively. Putting these values into (45) gives a ~ 1CT 10 which is 
in good agreement with the value obtained by Milgrom. It seems similar calculation for other 
galaxies gives almost the same order of magnitude for the value of ao- 
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8 Conclusion 

The problem of rotation curves for the galaxies has been investigated in the context of a 
new model which may be considered as an alternative to the dark matter and MOND. This 
model is motivated by the so called Weizsacker's Liquid- Droplet model which had been used to 
describe the structure of the nucleus, according to which the nuclear force is saturated inside 
the nucleus to exhibit a short range force and a flat curve of binding energy. Similarly, we 
have tried to show that the reason for the flat rotation curve of a galaxy may have its origin 
in the fact that the gravitational force exhibits an effective local short range interaction in the 
galactic structure while representing characteristic features of Mach principle. This kind of 
gravitational interaction at the galactic scale gives the rotation curve in good agreement with 
observations. We have also developed the Virial theorem and also obtained the Tully-Fisher 
relation. Moreover, we have given a physical explanation for the so called constant acceleration 
in the MOND. 
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